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Aiming to remedy the incorrect asymptotic behavior of conventional semilocal exchange- 
correlation (XC) density functionals for finite systems, we propose an asymptotic correction (AC) 
scheme, wherein an exchange density functional whose functional derivative has the correct (— 1/r) 
asymptote can be directly added to any semilocal density functional. In contrast to semilocal density 
approximations, our resulting exchange kernel in reciprocal space exhibits the desirable singularity 
of the type 0( — l/q 2 ) as q — >• 0, which is a necessary feature for describing the excitonic effects in 
non-metallic solids. By applying this scheme to a popular semilocal density functional, PBE [J. P. 
Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)], the prediction of prop- 
erties that are sensitive to the asymptote, such as the highest occupied molecular orbital energies 
and Rydberg excitation energies, is significantly improved, while the prediction of properties that 
are insensitive to the asymptote, such as reaction energies and valence excitation energies, remains 
essentially the same as PBE. In addition, without loss of accuracy, two closely related AC schemes 
are developed for the efficient treatment of large systems. 



Over the past two decades, Kohn-Sham density func- 
tional theory (KS-DFT) Q has been one of the most pow- 
erful theoretical methods for the ground-state properties 
of large electronic systems. Its time-dependent extension, 
time-dependent density functional theory (TDDFT) 0- 
0] has gradually become popular for the study of excited- 
state and time-dependent properties. 

In KS-DFT, the exact exchange-correlation (XC) den- 
sity functional E xc [p] remains unknown and needs to be 
approximated. Accurate density functional approxima- 
tions to E xc [p] have been successively developed to ex- 
tend the applicability of KS-DFT to a wide range of 
systems. Despite the recent advances in the orbital- 
dependent density functional approach [f| , semilocal den- 
sity functionals remain popular due to their computa- 
tional efficiency for large systems and reasonable accu- 
racy for applications governed by short-range XC effects 
[f|. However, due to the associated several qualitative 
failures, semilocal functionals can produce erroneous re- 
sults in situations where the accurate treatment of non- 
locality of the XC hole is important • 

One of the important and long-standing subjects in 
KS-DFT is the asymptotic behavior of the XC potential 
Vxc(r) = SE xc (r)/5p(r). For finite systems, the exact 
v X c(r) is known to exhibit the correct Coulombic (—1/r) 



However, due to the severe self- 
14j , the XC potentials derived from con- 



decay as r — » oo 
interaction error 
vcntional semilocal functionals fail to describe the correct 
asymptotic behavior, yielding qualitatively incorrect pre- 
dictions for properties where the asymptotic region of the 
XC potential is crucially important, such as the vertical 
ionization potentials and high-lying (Rydberg) excitation 
energies of atoms and molecules [15|, [l6[ . 

Without resorting to the expensive orbital-dependent 
density functional approach, currently, perhaps the most 



popular approach in practice to improve the asymptote 
of the XC potential is the model XC potential approach 
17H25{. wherein an asymptotically corrected XC poten- 



tial is directly modeled. However, as a number of popular 
model potentials are found not to be functional deriva- 
tives, these model potentials are not variationally stable, 
and the associated XC energies are not uniquely defined. 
Therefore, properties obtained from these model poten- 
tials need to be carefully interpreted |26| . 

By contrast, for a system of N electrons, the Fermi- 
Amaldi (FA) XC functional [27j], 
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which is just ( — 1/iV) times the Hartree energy functional, 
appears to be the simplest XC functional whose func- 
tional derivative has the correct (—1/r) asymptote. How- 
ever, there are several problems with the use of the FA 
model [28(. While the FA XC potential is correct in the 
asymptotic region, it is inaccurate elsewhere. Besides, 
due to its delocalized XC hole, the FA model is not size- 
consistent, where the energy of a system composed of two 
or more well-separated subsystems is not identical to the 
sum of the energies of the separate subsystems [2^] . 

In this Letter, we propose an asymptotic correction 
(AC) scheme for any system composed of atoms (e.g., 
atoms, molecules, and solids), wherein a modified FA 
XC functional, which is size-consistent in the calculated 
energy and whose functional derivative has the correct 
(—1/r) asymptote, can be directly added to any semilocal 
functional. Without loss of accuracy, two related efficient 
schemes are also developed for large systems. 

By partitioning and localizing a modified FA XC hole 
to the atoms in a system, we propose the "localized" FA 
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(LFA) exchange functional, 
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to resolve the size-inconsistency issue associated with the 
FA model [3(| . Here the second sum is over all the atoms 
in the system, PaA t ) is the cr-spin (a — a for spin up 
or j3 for spin down) electron density associated with the 
atom A, 



PaA v ) = WA(r)p a (r), 



(3) 



and the weight function waM, ranging between and 1, 
is of the Hirshfeld type 3jJ, |32J : 



HbPW 
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where p\{v) is the spherically averaged electron density 
computed for the isolated atom A. Na,<j is the number 
of the cr-spin electrons associated with the atom A, 



Na,<j = / pA,a{r)dr, 



(5) 



and the long-range interelcctron repulsion operator 
erf (a; |r — r'|)/ |r — r'| is to retain the correct asymptotic 
behavior without the (unneeded) energy contribution 
from the complementary short-range operator, where to 
is a parameter defining the range of the operators. Due 
to the sum rule of J^a w a{y) = 1, we have J2a PaA t ) = 
J2 A WA{r)p a {r) = p a (r) and Y,a N a,° = n <j ( tne num- 
ber of cr-spin electrons). 

By taking the functional derivative of E^ FA [p a ,pp], the 
LFA exchange potential for cr-spin electrons is 
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In the asymptotic limit, v^ FA (r) has the correct asymp- 
totic form, 
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From Eqs. and ([5]), we have 

E^ A [ Pa ,p f3 ] = 5 E / PAr)vZ A (r)dr, 



(8) 



which shows that the LFA exchange energy density per 
electron also has the correct (— l/2r) asymptote [17tl33j|. 

For calculations of excitation energies from adiabatic 
linear- response TDDFT H, the functional derivative of 
v x CT A ( r ) is taken to give the LFA exchange kernel for cr- 
spin electrons, 
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In contrast to semilocal density approximations, the LFA 
exchange kernel in reciprocal space has the correct long- 
wavelength 0{— l/q 2 ) divergence as q — > 0, which is cru- 
cially important for the pro per description of excitonic 
effects in non- metallic solids [34hi37|. We emphasize that 
this striking feature appears naturally from our fully non- 
local LFA exchange functional. 

To improve its description of short-range XC effects, 
the LFA exchange functional is combined with a popu- 
lar semilocal functional, PBE [38j | . However, this will 
produce a double-counting (DC) energy £l)c, which 
needs to be removed. As a significant fraction of the 
LFA exchange energy should be from the core regions 
of atoms, we presuppose that pA,a( r ) is strictly localized 
at Ra (the position of atom A) in Eq. @, which gives 
PaA v ) ~ ^A,a5( r — R-a) for satisfying Eq. ([5]), to esti- 
mate -EdCj 
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This estimate is very accurate for systems with highly 
localized charges (e.g., HC1), and less accurate for sys- 
tems with delocalized charges (e.g., benzene) [30]. Our 
resulting LFA-PBE functional is given by 
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Note that v^(r) = SEnc/Spai?) = —w/y/n is simply a 
constant. As the KS potential is only defined within an 
arbitrary constant, without loss of generality, we require 
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the KS potential to vanish asymptotically, which gives 
D°°(r) = 0. Unlike the FA model, LFA-PBE is size- 
consistent for any system composed of atoms 30]. Note 
that LFA-PBE (with ui = 0) reduces to PBE. 

For systems composed of many atoms, LFA-PBE can 
be computationally unfavorable due to the numerical 
integration of many Hartree-like potentials in v^ F J^(r). 
To resolve this computational bottleneck without loss of 
much accuracy, ^ FA (r) can be efficiently evaluated based 
on the resolution-of-identity (RI) approximation 3^, 4(| . 
Following Ref. 40], pA.a{ Y ) is first expanded with an aux- 
iliary basis set {g p (r)}, 



PaA v ) « p A A v ) = ^2a p g p (r), 



(12) 



where the expansion coefficients {a p } are given by Eq. (8) 
of Ref. [!o| (with p(r) being replaced by p A A r ))- The 
RILFA exchange potential is evaluated by 



RILFA 



(r) 



E 



N A J PAAr) |r-r'| 



(13) 

From Eq. (10) of Ref. [40j (with the Hartree energy be- 
ing replaced by the LEA exchange energy), the RILFA 
exchange energy is given by 
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The asymptote of Ae (r) is easily shown to be correct 
hm ^A s(r) = _£ WA(r) 1 

A ' ' 

= --E WA(r ) = 

r r 
A 

Based on Eq. ©, the LFAs exchange energy is given by 
^ LFAS = \ E / Ar(r)< FAs (r)dr, (17) 

a=a,f3 

to retain the correct (— l/2r) asymptote of the LFAs ex- 
change energy density per electron. Although w^ As (r) 
differs from the functional derivative of E^ FAs by a factor 
of 2, the prescribed LFAs scheme approaches to the LEA 
scheme for a sufficiently small uj value, where v^ FAs (r) 
becomes an excellent approximation of v^ FA (r). Simi- 
larly, LFAs-PBE is defined by Eq. (HT), with E% FA being 
replaced by ^ FAs . 

For the exact KS-DFT, the minus HOMO energy of 
a molecule should be the same as the vertical ionization 
potential (IP) of the molecule 0, [H EH-S] . Therefore, 
the optimal uj values for LFA-PBE, RILFA-PBE, and 
LFAs-PBE are determined by fitting the predicted IPs 
(calculated by the minus HOMO energies) of 18 atoms 
and 113 molecules in the IP 131 database to the corre- 
sponding experimental IPs [44[ . All calculations are per- 
formed with a development version of Q-Chem 3.2 [451 ]. 
using the 6-311++G(3df,3pd) basis set (and sufficiently 
large auxiliary basis sets for the RILFA scheme), unless 
noted otherwise. The error for each entry is defined as 
(error = theoretical value — reference value). 



For a sufficiently large auxiliary basis set, the RILFA 
scheme approaches to the LEA scheme. Here, we define 
RILFA-PBE by Eq. (TTTJ), with E% FA being replaced by 
£* ILFA . 

For very large systems, both the LEA and RILFA 
schemes may be impractical, compared to semilocal func- 
tionals. Aiming to retain the correct (— 1/r) asymptote 
with essentially no added computational cost relative to 
semilocal functionals, we define the LFAs exchange po- 
tential v^ FAs (r) ("s" refers to "strict"), by presupposing 
the strict localization of paA t ) at Ra in u^ FA (r) (i.e., 
PaA*) ~ N A ,*6{r-R A )), 

tiu(r) 
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FIG. 1. The root-mean-square (RMS) errors of the IP131 
database 0, for LFA-PBE, RILFA-PBE, and LFAs-PBE at 
different values of to. The to = case corresponds to PBE. 
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Experimental ionization potential (eV) 



TABLE I. Mean absolute errors (in kcal/mol) of the reaction 
energies of 30 chemical reactions [4a ]. calculated by the PBE 
and LFA-corrected PBE functionals (uj = 0.15 Bohr -1 ). (1 
kcal/mol = 0.0434 eV = 0.00159 hartree.) 



System 


PBE 


LFA-PBE 


RILFA-PBE 


LFAs-PBE 


Reaction (30) 


4.38 


4.48 


4.47 


4.42 



TABLE II. Mean absolute errors (in eV) of the 19 valence and 
23 Rydberg excitation energies of 5 molecules [4^], calculated 
by the PBE and LFA-corrected PBE functionals (w = 0.15 
Bohr- 1 ). 



System 


PBE 


LFA-PBE 


RILFA-PBE 


LFAs-PBE 


Valence (19) 


0.32 


0.29 


0.29 


0.29 


Rydberg (23) 


1.30 


0.46 


0.46 


0.49 



FIG. 2. Calculated versus experimental ionization potentials 
(IPs) for the IP131 database Q. Inset I shows the differences 
between the IPs calculated by RILFA-PBE and LFA-PBE, 
while inset II shows the differences between the IPs calculated 
by LFAs-PBE and LFA-PBE. For the LFA schemes, u = 0.15 
Bohr" 1 is adopted. 



As shown in Fig. [TJ the minimum root-mean-square 
(RMS) errors of the IP131 database, for LFA-PBE, 
RILFA-PBE, and LFAs-PBE all occur at u = 0.15 
Bohr -1 , where the RMS errors associated with the LFA- 
corrected PBE functionals are more than three times 
smaller than the RMS error associated with PBE (the 
uj = case), reflecting the importance of the correct 
asymptote of XC potential for the accurate prediction 
of IPs 0. 

Adopting uj = 0.15 Bohr -1 for the LFA schemes, the 
calculated IPs are plotted against the experimental val- 
ues in Fig. [2j As can be seen, the differences between the 
IPs calculated by RILFA-PBE and LFA-PBE are within 
0.005 eV, and the differences between the IPs calculated 
by LFAs-PBE and LFA-PBE are within 0.1 eV. There- 
fore, the LFA-corrected PBE functionals yield very simi- 
lar results, which implies that their XC potentials should 
be very similar [30]. By contrast, the predicted IPs of 
PBE are seriously underestimated, due to the incorrect 
asymptote. 

Similar results are found, when our LFA schemes are 

combined with LDA (UBS]. As both the LDA and PBE 
XC potentials decay exponentially in the asymptotic re- 
gion, their predicted IPs remain similar, requiring essen- 
tially the same corrections from the LFA schemes (i.e., 
with the same optimal uj) [30]. Therefore, uj — 0.15 
Bohr -1 can be recommended as the optimal uj value, 
when the LFA schemes are combined with a local or 
semilocal functional whose functional derivative has the 
(incorrect) exponential asymptote. 

The performance of the LFA schemes is also examined 



on the reaction energies of 30 chemical reactions (a test 
set described in Ref. [3]). As shown in Table U the 
performance of the LFA-corrected PBE functionals (uj = 
0.15 Bohr -1 ) is very similar to that of PBE do|. As 
reaction energies are rather insensitive to the asymptote 
of the XC potential, our schemes do not affect the good 
performance of PBE. 

To assess the performance of our LFA schemes on va- 
lence and Rydberg excitation energies, we perform adi- 
abatic linear-response TDDFT calculations, using the 6- 
311(2+, 2+)G** basis set, on five molecules: nitrogen 
gas (N 2 ), carbon monoxide (CO), water (H 2 0), ethylene 
(C2H4), and formaldehyde (CH 2 0) on the experimental 
geometries taken from Ref. [49j . For computational sim- 
plicity, we adopt /^ A (r,r'j = (rather than Eq. ©) 
in our TDDFT calculations. For finite systems, this ap- 
proximation should not make much difference in the pre- 
dicted excitation energies 17j. As shown in Table HH the 
LFA-corrected PBE functionals (a; = 0.15 Bohr -1 ) per- 
form well for both the valence and Rydberg excitations, 
while PBE severely underestimates Rydberg excitation 
energies due to the lack of the correct asymptote of the 
XC potential 

In conclusion, we have developed the LFA scheme, 
wherein an exchange density functional whose functional 
derivative has the correct (— 1/r) asymptote can be di- 
rectly added to any semilocal density functional. In con- 
trast to semilocal density approximations, the LFA ex- 
change kernel in reciprocal space is shown to exhibit the 
desirable singularity of the type 0(—l/q 2 ), which is an 
important feature for the description of excitonic effects 
in non-metallic solids. Applying the LFA scheme to PBE, 
the resulting LFA-PBE {uj = 0.15 Bohr -1 ) has yielded 
accurate IPs and Rydberg excitation energies for a wide 
range of atoms and molecules, while performing similarly 
to PBE for properties insensitive to the asymptote, such 
as reaction energies and valence excitation energies. In 
addition, without loss of accuracy, two closely related 
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schemes (RILFA and LFAs) have been developed for the 
efficient treatment of large systems. As the predicted 
results of the LFA-corrected PBE functionals are very 
similar, LFAs-PBE (lu = 0.15 Bohr -1 ), which has the 
correct (— 1/r) asymptote with essentially no added com- 
putational cost relative to PBE, is potentially very useful 
for the study of the frontier orbital energies, excitation 
energies, and time-dependent properties of large systems. 
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